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Abstract 

We consider the Cauchy problem for the massive Dirac equation 
in the non-extreme Kerr-Newman geometry outside the event horizon. 
We derive an integral representation for the Dirac propagator involving 
the solutions of the ODEs which arise in Chandrasekhar's separation 
of variables. It is proved that for initial data in L"^^ near the event 
horizon with decay at infinity, the probability of the Dirac particle 
to be in any compact region of space tends to zero as t goes to infinity. 
This means that the Dirac particle must either disappear in the black 
hole or escape to infinity. 
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1 Introduction 



It has recently been shown that the Dirac equation does not admit normaliz- 
able, time-periodic solutions in the non-extreme Kerr-Newman axisymmetric 
black hole geometry [1]. This was interpreted as an indication that a Dirac 
particle either falls into the black hole or escapes to infinity, but that it can- 
not stay in a bounded region outside the event horizon. In this paper we 
make precise this interpretation in the general time-dependent setting. We 
thus consider the Cauchy problem for the Dirac equation with smooth ini- 
tial data on the hypersurfacc t = 0, compactly supported outside the event 
horizon. We study the probability for the Dirac particle to be inside a given 
annulus located outside the event horizon, and we prove that this probability 
tends to zero as t goes to infinity. Hence, in contrast to the situation for a 
bounded orbit of a classical point particle, there exists no compact region 
outside the event horizon in which the quantum mechanical Dirac particle 
will remain for all time. In more precise form, our result is stated as follows. 
Recall that in Boyer-Lindquist coordinates {t,r,'d,(p) with r>0,0<t?<7r, 
< If < 27r, the Kerr-Newman metric takes the form [2] 

ds^ = gjk dx^x^ 

= §{dt - a sin^ ^ - u(^^+ (1.1) 
- ^ (a dt - (r^ + a^) d^f 

with 

U{r, d) = + cos^ ^ , A(r) = - 2Mr + + , 

and the electromagnetic potential is given by 

Q r 

Aj dx-' = — — {dt — a SIT? ■& dip) , 

where M, aM and Q denote the mass, the angular momentum and the 
charge of the black hole, respectively. Wc shall here restrict attention to 
the non-extreme case > -\- Q^. Then the function A has two distinct 

zeros, 

ro = M - - a2 - Q2 and n = M + M"^ - - Q'^ , 

corresponding to the Cauchy and the event horizon, respectively. We will 
here consider only the region r > r\ outside of the event horizon, and thus 
A > 0. 

Theorem 1.1. Consider the Cauchy problem for the Dirac equation in the 
non-extreme Kerr-Newman black hole geometry outside the event horizon 



{ijWj -m)'i>{t,x) = 0, ^-(0,^) = *o(a;) , (1-2) 
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where the initial data ^'o is in L^((ri, oo) x S*^, dfj,)^, where dfi is the induced 
invariant measure on the hypersurface t = const. Then for any 6 > and 
R > ri + S, the probability for the Dirac particle to be inside the annulus 
Kgn = {ri + 5 < r < R} tends to zero as t ^ oo, i.e. 

lim / {^-i^^){t,x)v. djjL = 0, (1.3) 

where v denotes the future directed normal. 

The decay of probabilities in compact sets (1.3) can be stated equivalently 
that the Dirac wave function decays in Lf^^ outside and away from the event 
horizon. Since the Dirac equation is linear and stationary, for smooth initial 
data we obtain immediately that also the time-derivatives decay in 

L^p^,, and standard Sobolev methods yield that ^ decays even in L^^. We 
point out that the initial data is merely bounded (but not necessarily small) , 
near the event horizon. Our assumptions include the case when the initial 
data is smooth and bounded in the maximal Kruskal extension up to the 
bifurcation 2-sphere (as is considered in [7] for the wave function in the 
Schwarzschild geometry). We note that the axisymmetric character of the 
background geometry makes the analysis significantly more delicate than in 
the spherically symmetric case, mainly because for a 7^ both the radial and 
angular ODEs depend on the energy, and thus for the study of the dynamics 
we must consider the system of these coupled equations. 

The proof is organized as follows. We first bring the Dirac equation 
into the Hamiltonian form idt^ = with a self-adjoint operator H. Our 
main technical tool is an integral representation for the Dirac propagator 
exp(—itH) acting on wave functions with compact support outside the event 
horizon. This integral representation is stated in Theorem 3.6. To derive it, 
we first consider the Dirac equation in an annulus outside the event horizon 
with suitable Dirichlet-type boundary conditions. The Hamiltonian corre- 
sponding to this modified problem has a purely discrete spectrum, and thus 
the propagator can be decomposed into discrete eigenstates. We then take 
the infinite- volume limit by letting the inner boundary of the annulus tend 
to the event horizon and the outer boundary to infinity in a suitable way. 
Our construction is based on Chandrasekhar's separation of variables for the 
Dirac equation in the Kerr-Newman metric [3, 4, 5] together with estimates 
for the asymptotic behavior of the amplitudes and phases of the separated ra- 
dial cigcnfunctions (Lemmas 3.1 and 3.5), and for the spectral gaps (Lemma 
3.3). The usefulness of our integral representation for the propagator is that 
it explicitly gives the continuous spectral measure of H in terms of the solu- 
tions of the radial and angular ODEs arising in Chandrasekhar's separation 
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of variables. For initial data which is compactly supported outside the event 
horizon, the decay of the probabilities (1.3) then follows by standard results 
of Fourier analysis. The generalization to initial data in and near 
the event horizon is done by approximation in our Hilbert space framework. 



2 Separation of Variables, Hamiltonian Formula- 
tion 



The Dirac equation in the Kerr-Newman geometry can be completely sep- 
arated into ODEs by Chandrasekhar's method [3, 4, 5]. We here outline 
the separation, see [1] for details. After the regular and time-independent 
transformation 

^ ^ 4f = 



with 

S = A4 diag ^(r — za cosi?)2, (r — za cos??)2, 

(r + ia cos t?) 2 , (r + ia cos i})^^ , 

the Dirac equation can be written as 

+ = 



with 



7^ 



A 



imr 










\ 







-imr 
















—imr 









V 







imr 


/ 




—am cost? 

















am cos?? - 


£_ 












—am cos?? 

















am cos 



(2.1) 



(2.2) 



and the differential operators 

P± - — T - 

dr A 



(r^ + a^)!- + a-- 



dt 



d_ 

dip 



ieQr 



d cot I? 



' . a 9 Id 
a sm??— + — — — 
ot smw Oip 
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Employing the ansatz 



/ X_{r)Y_{'d) \ 

X+{r) Y-i'd) 
V X_(r)y+(t?) / 



keZ, (2.3) 



we obtain the eigenvalue problems, 
7^* = A#, 



-A* , 



(2.4) 



under which the Dirac equation (2.2) decouples into the system of ODEs 

= (2.5) 
= 0, 



\/A X)+ imr - A \ / X+ 
-imr-X yfKV_ ) \ X_ 

£+ — am cos ■!? + A 

am cos 1? + A —C- 



y+ 



where V± and C± are the radial and angular operators 



d i 



u (r^ + 



k + -] a + eQr 



d cot 7? 



au smv + 

smv 



(2.6) 

(2.7) 
(2.8) 



We will in what follows also use the vector notation X = (X+,X_), Y = 
(y_,y+) and for clarity sometimes add indices for the parameters involved, 
e.g. X'^^^ = X. We point out that (2.3) is an eigenfunction of the angular 
operator id^p with eigenvalue k + ^. The reason why we need to consider 
half odd integer eigenvalues is that the transformation from the usual single- 
valued wave function in space-time, to the wave function ^ in (2.2) involves 
a sign flip at = (see [1, Section 2.1]). 

In this paper, wc want to study timc-dcpcndcnt solutions of the Dirac 
equation. In the separation ansatz (2.3), the dynamics of the solution is 
encoded through the w-dependcncc in the ODEs (2.5) and (2.6). Unfortu- 
nately, both the radial and angular operators (2.7) and (2.8) depend on w, 
making the situation rather complicated. Therefore it is useful to bring the 
Dirac equation (2.2) into Hamiltonian form, in a way which is compatible 
with the separation of variables. We first bring the time derivative in (2.2) 
to one side of the equation and obtain 

B + asiu^C] i — ^ = {V? + A^)^ (2.9) 



A 
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with 



B = 



/ -i 0\ 

i 

1 

V -i y 



c = 



/ 



-1 

V -1 








-1 \ 




/ 



where the operators TZ'^ and A'^ are obtained from TZ and A by setting 
the time derivatives to zero. The matrices B and C satisfy the relations 
B^ = 1 = and BC = CB. Thus the hnear combination of these matrices 
which appears in (2.9) can be inverted with the formula {aB + (3C)~^ = 
{a^ - (aB - (5C) (and a,/3 e IR). Furthermore, we introduce a new 

radial variable u€ (—00,00) by 



du 
dr 



+ 0? 



(2.10) 



Omitting for simplicity the hat of the wave function, the Dirac equation (2.(; 
becomes 



• 5 T 



with the Hamiltonian 
■(r2 + a2)2 



H = 



A 



- sin^ ■!? 



"1 / 2 I 2 



1 - 



A sin^ i9 
(r^ + a 



2^2 



1 



A 

a \/A sin 'd 



B - a siwdC] {n^ + A^) 



+ 



BC 



(7^ + ^), (2.11) 



where r is an implicit function of u, and 

mr \/A 



7^ = - 



^2 _j_ q2 



+ 



A 



am cos i9 V A 

j,2 _|_ q2 










/ 1 \ 
1 
10 

V 1 / 

/ i \ 

i 

-i 

V -z / 

-X+ \ 



M 

M_ / 



^_ 














£+ 














£+ 














-£ 



(2.12) 



(2.13) 
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with 



d { ia d eQr 



du yr"^ + a? d(p + a 

y/A r d cot 1? Id 
M± = ^r- — n I i — + * 



+ \ 2 sin 1} dip 



The Hamiltonian (2.11) is an operator acting on the wave functions on the 
hypersurfaces t = const. The simplest scalar product on such a hypersurface 
is 

/oo rl r2n 

du dcosi) dip^it,u,'d,ip)^it,u,'d,ip) , (2.14) 
-00 J-i Jo 

where denotes the complex conjugated, transposed spinor. In the spher- 
ically symmetric case a = 0, the Hamiltonian (2.11) is Hcrmitian (i.e. for- 
mally self-adjoint) with respect to this scalar product. However for a ^ 0, 
H is not Hermitian. In order to get around this problem, we introduce 
a different scalar product as follows. Notice that the operators 1Z and A, 
(2. 12), (2. 13), are both Hermitian with respect to (2.14). The reason why 
the Hamiltonian (2.11) is not Hermitian is that, when the taking the adjoint 
of H using integration by parts, one gets r- and T^-derivatives of the square 
bracket in (2.11). But we can compensate this square bracket by inserting its 
inverse into the scalar product. Thus we introduce on the four-component 
spinors the inner product 



<*l^>(t,M,^) = ^it,u,^,ip) (1 + ^^^-^Bc]^t,u,^,ip) (2.15) 



and define the scalar product <.|.> by 



du dcosT? / dip \^>(^t,u,^,cp) ■ (2.16) 
-00 J-i Jo 

Then the Hamiltonian H is Hermitian with respect to (2.16). Let us verify 
that (2.16) is positive. In the region outside the event horizon under consid- 
eration, r > ri > M. Also, since we are in the non-extreme case, M > Q,a, 
and as a consequence, A < r^. We conclude that 



a s/A sin 



J.2 _|_ q2 



a^/A a ^/A 
_j_ r r 



Combining this inequality with the fact that the matrix BC has eigenvalues 
±1, we obtain that the bracket in (2.15) is indeed a positive matrix. 
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Wc denote the Hilbert space of wave functions with scalar product (2.16) 
by H. Then the operator H, (2.11), is essentially self-adjoint on H with 
domain of definition 

D{H) = C^(]R X . 

In Section 3, we shall consider the Dirac operator also with certain Dirichlet- 
type boundary conditions, which we now introduce. First for given U2 G 
IR, we restrict u to the half line u € (—00,^2] and impose the boundary 
conditions 

^'i(u2,T?,¥') = 5'3(M2,t?,¥5) and ^'2(ii2, 1?, '^) = *4(^i2, 1?, '^) • (2.17) 

Let 7Y„2 be the Hilbert space of square integrable wave functions '^{u, ip), 
u < U2 with the scalar product 

du dcos-d dip <^\^>^t,u,d,v) ■ (2-18) 
-00 J-i Jo 

Then the Hamiltonian (2.11) on TC^^ with boundary conditions (2.17), which 
we denote for clarity by Hu2, is Hermitian (the main point here is that the 
boundary values at u = U2 vanish when the adjoint of Hy,^ is calculated 
using integration by parts). The operator Hy,^ is essentially self-adjoint on 
Hu2 with domain of definition 

D{Hu^) = {-^e C^((-oo,U2] X 5^)^ and (2.17) is satisfied} . 

Similarly, for ui,U2 G IR, ui < U2, we restrict u to the closed interval 
u G [^1,^2] with boundary conditions 

*i(«i) = *3(wi), *2(wi) = *4(wi) ^2 19) 

*l(«2) = *3(W2), ^2(^2) = ^4(^2) . ^ ' ' 

Wc denote the Hilbert space of square integrable wave functions ^'(u, 1?, (^), 
ui ^ u ^ U2, with the scalar product 

<^>ui,u2 ■■= / du dcosi? / dip <'i>\^>^t,u,i},v) (2-20) 

Jui J-1 Jo 

by Ti.ui,u2J and the Hamiltonian (2.11) on 'Hui,u2 by -ffui,u2- It is essentially 
self-adjoint with 

-D(-ff„i,„2) = G C^{[ui,U2] X S^f and (2.19) is satisfied} . 

Our above Hamiltonian formulation of the Dirac equation is well-suited 
to Chandrasekhar's separation of variables. Namely, the boundary condi- 
tions (2.17) reduce to simple boundary conditions for the radial functions. 



X+{U2) = X_{U2) , 



(2.21) 
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whereas (2.19) amounts to 

X+{ui) = X_{ui) and X+{u2) = X_{u2) . (2.22) 

The scalar product (2.14) sphts into the product of a radial and an angular 
part, namely 

with 

/oo 
Xk^^{u) X''''^'^'{u) du 
-oo 

The scalar product (2.16), however, does not split into a product, more 
precisely 



+a (X'^'^^ I I x'^''^'^') (y*^'^^ I sint? £7^ I y'^''^'^') , (2.23) 

J.Z _|_ 

where cr* are the Pauli matrices. This mixing of the radial and angular parts 
in the scalar product can be understood from the fact that the Kerr-Newman 
solution is only axisymmetric. 



3 An Integral Representation for the Propagator 

The propagator exp{—itH) has the spectral decomposition 

/oo 
e-'^^ dE^ , (3.1) 
-oo 

where dE^^ is the spectral measure of H. In this section, we shall bring 
this formula into a more explicit form. This will be done by expressing 
the spectral measure in terms of solutions of the radial and angular ODEs 
of the previous section. Since the spectrum of H is continuous, it is not 
obvious how to relate the spectral measure to the solutions of our ODEs. 
To bypass this problem, we begin with the spectral decomposition of the 
operator Hui,u2 (which has a purely discrete spectrum), and then deduce 
the desired integral representation for exp(— ziif) by taking suitable limits 
ui — — oo and U2 — oo. 
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As an elliptic operator on a bounded domain, the Hamiltonian Hui,u2 
has a purely discrete spectrum with finite-dimensional eigenspaces (see [6]). 
In view of our separation of variables, the most convenient eigenvector basis 
is the following. First we can choose the basis vectors as eigenvectors of the 
operator id^p with eigenvalue k + ^, k ^ Z. We denote this eigenspace of idip 
by H^^ u^, and the restriction of Hui,u2 to Ti-^^ ^^ by H^_^^^^. Furthermore, 
the basis vectors can be chosen as eigenvectors of the angular operator A. 
As is shown in the Appendix, the spectrum of A on "H^j „2 discrete, non- 
degenerate, and depends smoothly on lu. Thus the eigenvalues of A can be 
written as A„(cl'), n € Z, with A„ < Xn+i and A„(.) G C°°(IR,). For any given 
k & Z, io G c(-f^ui,u2)' n G Z, the radial ODE (2.5) has at most one 
solution satisfying the boundary conditions (2.19). Hence we have for any k, 
u, and n at most one eigenstate of H^-^^^u^, which we denote by The 
set of n for which such an eigenvector exists is denoted by N{k,Lo). Thus 
our eigenvector basis is 

We normalize these eigenfunctions with respect to the scalar product (2.14); 
more precisely, we normalize both the radial and angular parts according to 

{xtru2 1 xtru.) = 1 ' (^'"" I = i (3.3) 

with X and Y as in (2.3). Since the angular operator A is self-adjoint with 
respect to the scalar product (.|.), its eigenvectors are orthogonal, and thus 
the eigenfunctions for fixed k and lo are even orthonormal, 

(*Sr«2 I <Tu.) = ' e N{k,u;). (3.4) 

We mention for clarity that for different values of uj, the eigenfunctions are in 
general not orthogonal with respect to (.|.), but since i^ui,u2 is self-adjoint 
with respect to <.|.>, its eigenspaces are orthogonal with respect to the 
latter scalar product, and thus 

These subtle differences between the two scalar products clearly become 
irrelevant in the spherically symmetric case a = 0. 

In the basis (3.2), the spectral decomposition (3.1) for -ffui,u2 can be 
written as 

= E E E w<:2<<:2i^>)-(3-5) 

keZ uj&a{Hk^^^^) \n,n'&N{k,uj) ) 



(3.2) 
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Here the coefficients c„„' must be chosen such that the bracket in (3.5) is the 

jk 



projection of ^' onto the eigenspace of Hj^ corresponding to the eigenvalue 



u; more precisely, 

Cnn' = {A-X^, With ^nn' = <*t«2l*™>- M 

Notice that the first two sums in (3.5) give a decomposition of into the 
orthogonal eigenstates of the operators ^^^p and H, respectively, and thus 
converge in norm in 'Hui,u2- The bracket in (3.5) is the basis representation 
of the projector on the respective eigenspace. 

Our first goal is to take the limit ui —oo in (3.5). We expect that 
in this infinite volume limit, the "energy gaps" ALOkn between neighboring 
eigenvalues, defined by 



Aukn = min{a)fe„ - LOkn I f^fcrt > ^^kn} with 

jk 



^kn,'^kn e (y{H^ ) and N{k,LUkn), N{k,iJkn) 7^ , 



should tend to zero. The basic idea is to rewrite the sum over the spectrum 
in (3.5) as Ricmann sums which converge to integrals as ui —oo, yielding 
a formula for the propagator of the Hamiltonian Hu2- For making this idea 
mathematically precise, it is essential to get good estimates for Au^n aiid to 
relate the eigenvectors ^u^^2 (^•^) solutions 

W with k ez, u eM, nez, ue (-00, U2] 

of the Dirac equation with boundary conditions (2.17). We note that the 
convergence of the series in n is not a real issue. Indeed, using an ap- 
proximation argument, we will show in the proof of Theorem 3.6 that we 
may restrict attention to a finite number of angular momentum modes. 

We denote the radial and angular functions corresponding to ^'fj^f" by 
X^^" and Y'^'^"-, respectively. In the variable u, (2.10), the radial equation 
(2.5) becomes 



d^+^'^^^H -1 



X = ^ ( ° imr - A , ^ 
+ V —imr — X ' 



with 



, {k + ^)a + eQr 



and where for ease in notation the indices of X were omitted. The next 
lemma describes the asymptotic behavior oi X{u) as n — — oo. 
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Lemma 3.1. Every nontrivial solution X of (3.7) with boundary conditions 
(2.21) is asymptotically as u ^ — oo of the form 



Xiu) = ( ^,n,u}^ ] + Mu) (3.8) 



with 



fo + (3.9) 

(fc + ^)a + eQri 

= Lo ^ ^ — (3.10) 

r{ + 

l-Rol < ce''" (3.11) 
and suitable constants c,d>Q, which can be chosen locally uniformly in to. 



Proof. Substituting into (3.7) the ansatz 

-iUou 



Xiu) 



we obtain for / the equation 



(3.12) 



du 



A 



1 
-1 



+ 



j,2 _j_ Q,2 



e-2^^"(imr - A) 

,2inu^_-^^ - A) 



/. (3.13) 



The square bracket vanishes on the event horizon r = ri. In the variable 
?i, this leads to exponential decay for n — ^ — oo, in the sense that there are 
constants ci , d > such that 



du 



< ci e'^" I/I 



(3.14) 



Since X is a nontrivial solution, |/| ^ 0. Thus we can divide (3.14) by |/| 
and integrate from any u < U2 to U2 to obtain 



log I/I 



< C2 e 



du 



U2 

u 



with C2 = ci/d. Since the right side of this inequality stays finite when 
u ^ — oo, we conclude that there is a constant L > with 



I < I/HI < L 



for all u < U2. 



(3.15) 
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Using that A depends smoothly on uj (see the Appendix), the constants 
ci, C2,d, and L clearly can be chosen locally uniformly in to. 



We substitute (3.15) into (3.14), 

< ciLe"'" . (3.16) 



du 



This inequality shows that /' is integrable, and thus f(u) converges for 
u — > — oo. Setting 

(3.15) 

/o = lim f{u) + , 

we can integrate (3.16) from — oo to -u < 6 and get 

\i{u)-h\ < ce"^ 

with c = c\Ljd. Substituting in the ansatz (3.8), we get (3.11). ■ 

^Prom (3.9) we see that X{u) does not decay to zero for u — oo. As a 
consequence, the function ^^"^ cannot have finite norm and thus is not a 
vector in the Hilbert space 'Hu2- This shows that the Hamiltonian has 
no point spectrum. In contrast to (3.3), we normalize the functions ^^j" 
according to 

lim \X^'f''\ = 1 , {Y^'''' I y*^'"") = 1 . (3.17) 

«— >— oo 

The next two lemmas describe the behavior of the normalization factors 
and the energy gaps as ui — > — oo. 

Lemma 3.2. For fixed U2 and asymptotically as ui — ^ — oo, 

xi^rru, = with (3.18) 

|5(ni)|-2 = iu2-ui) + 0(1) . (3.19) 



Furthermore, 



< <y'='^" I sin^a^ I y'"^'''> , (3.20) 

U2 - Ui 



where the constant c can be chosen locally uniformly in lo. 

Proof. Since X^'^J^^ and X^^"" are solutions of the same ODE (3.13) with the 
same boundary conditions at U2, they clearly coincide up to a normalization 
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factor g, i.e. Xjj^"^ ~ 9^U2^- Taking the norm on both sides and using the 
first part of (3.3), we obtain that 

We now substitute (3.8), multiply out, and use that |/op = 1 according to 
the first part of (3.17), to obtain 



ru2 

/ (1 + XRo + RqX - |i?oP) du . (3.21) 



Since X is bounded and Rq has exponential decay (3.11), the last three 
summands in (3.21) are integrable, and thus \g{u\)\~'^ — {u^ — Ui) is bounded 
uniformly in ui. This proves (3.19). 

The scalar product <^u'J',«2 I ^m%2> ^^"^ ^e computed via (2.23). The 
orthonormality (3.4) yields that 



= ci{X, 



kum 



a' I X^'^'^J (y'='^" I Sim? (7^ I y^^"') . (3.22) 



1 I \^k(jjn'\ 



In order to estimate the radial scalar product, we first note that the factor 
\/A goes exponentially to zero for u — oo, and thus 



kuin 



I vkui', 



2 I "v^fcwn' \ ^ „ I pdu I vkiAjn i i -y-kum' i j 



for some constant C4 > 0. Substituting (3.19) and using that the integral is 
uniformly bounded due to the factor exp du, we obtain the estimate 



kum 

y-i,y-2 I „2 



_2 I \rkujn' ^ 
\ ^Ul,U2> 



< C5 {U2 - ttl) ^ , 



which together with (3.22) yields (3.20). 



Lemma 3.3. The following estimate holds asymptotically as ui — —00, 

Awjkn = — - — + 0{{u2 - ui)-^) , (3.23) 

U2 — Ui 

for fixed U2 locally uniformly in to. 
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Proof. We consider solutions of (3.7) satisfying the boundary conditions 
at U2 and ask for which values of lo and A„(ix') our boundary conditions are 
also fulfilled at ui. As is immediately verified from (3.7), 

A(|X,|^_|X-P)=0. 

Thus {X^l"^ — is independent of u, and since it vanishes at U2, 

= forall«<U2. (3.24) 

Hence for the boundary values at ui, we need not be concerned about the 
absolute values of X± ; it suffices to consider the condition for the phases 

argX+(ui) = argX_(iti) . (3.25) 



It is convenient to work again with the ansatz (3.12). In order to describe 
the dependence on lo, we differentiate (3.13) with respect to lo. Since A 
depends smoothly on w, we obtain the bound 



d_ 
du 



du^f 



< ci e'^" \d^f\ + C3 e'*" I/I 



du 



(3.26) 



with constants ci, d as in (3.14) and C3 > 0. Using that |/| is bounded from 
above (3.15), we get 



< <!ie''"(|S»/l+<:4) 



with C4 = cj,Ljc\. Similar to the development after (3.14), dividing by 
(|9u;/| + C4) and integrating yields 



tl2 

u 



iog(|5^/i+c4)i::^ < C2e'^" 

and since the right side of this inequality is uniformly bounded in 

\d^j\ < C5 (3.27) 

for some constant C5 > 0. 

For the study of the phase shifts, we introduce the phase function 

p{u) = arg/+(n) - arg/_(u) - 2il ^2 

(the last summand was included so that p{u2) = 0). The derivative of the 
argument of a complex-valued function h is given by 
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Using this formula together with the fact that, according to (3.15) and (3.24) 
both |/_|_| and |/_| are bounded away from zero, we obtain that 



d_ 

du 



dujp 



< 4L 



d_ 

du 



du^f 



+ \d^f\ 



d_ 

du 



f 



Substituting in (3.26), (3.16) as well as the bounds (3.15) and (3.27), we 
conclude that 



du 



< C6 e 



du 



with some constant cg > 0. We integrate this inequality from u < U2 to U2- 
Since p{u2) = independently of lo, the boundary term dojp{u2) drops out, 
and we obtain the bound 



|5c.p(n)| < C 



for all u < U2 



(3.28) 



with a constant C > 0. This means that the equation for /, (3.13) leads 
only to finite phase shifts. 

The boundary conditions at ui, (3.25), are fulfilled iff 

$ := 20 {u2 - ui) + p = (mod 27r) . 

Differentiating with respect to u and integrating again from ujj to uju, ooj < 
uji, we obtain that 



\^{ujii)-^{(jji) - 2{ujn - u)i) {u2 - ui) 
< 



i^ii (3.28) 

\dujp\ du < C {ioii - ui) , 

and this proves (3.23). ■ 

Wc can now prove the integral representation for the propagator of ■ 
Proposition 3.4. For every ^ G Cq°((— oo, ■U2]) x S'^)^ and x = (u, (/p), 

1 /-oo 



Proof. According to the bound (3.20), the operator A in (3.6) converges 
uniformly in uj and k to the identity as ui — > —00, and thus (3.5) simplifies 
asymptotically to 



(e '*-'^"i.''2 x]/) (x) 

= E E ' 



-iu)t 



E^kujn 

neN{k,uj) 



<^™l*> + 0{{U2-U, 
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Using (3.18) and (3.19), we can express '^'^^^ by 



(e---.*)(x) = E7Z^ E 



Uo — U\ 

neN{k,u>) 

The gap estimate, Lemma 3.3, shows that the sum over the spectrum is a 
Riemann sum which converges as -ui — ^ — oo to an integral. ■ 



du 





—iuj 


im ^ 




[( 










—im 







X 



The idea for proving an integral representation for exp{—itH) is to take 
in (3.29) a suitable limit U2 +oo. In preparation, we need to derive 
estimates which describe the asymptotics of solutions of the radial equation 
(3.7) for large u. In this regime, 

—ieQ —iniM — A 
imM — A ieQ 

+ 0{u-^)X. (3.30) 

Thus the matrix potential on the right converges for u — > (X). If |w| < m, 
its eigenvalues A = ±Vm?^^uj^ are real, and this leads to one fundamental 
solution of (3.30) which decays exponentially like cxp(— -y/m^ — lj'^u), and the 
other solution has exponential growth ~ ex.p{\/'m'^ — uo'^u). We denote these 
two fundamental solutions by ^i^"^ and respectively, and normalize 

them according to 

lim \^\T{u)\ = 1 , (3.31) 

M— > — OO ' 

where our notation 1/2 means that the above equation is valid in both cases 
1 and 2. For |a;| > m, on the other hand, the eigenvalues of the matrix 
potential at n = oo are imaginary, A = ±i\/u>'^ — m^, and this leads to two 
fundamental solutions with oscillatory behavior ~ exp(zbi\/<^^ — w?u). 
For the normalization, we are now free to choose both the amplitude and 
the phase. Our convention is that 

/o'T = [l] /o^J = ( J ) (3.32) 

with /q'^ as in the asymptotic expansion (3.8). The next lemma describes 
the asymptotics of the oscillatory solutions as u ^ oo. 

Lemma 3.5. Every nontrivial solution X of (3.7) for \u;\ > m has for large 

u the asymptotic form 

f£ ) + Roo{u) (3.33) 
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with 



foe ^ (3.34) 

f o ^(Q + \ 

$ = e{uj) I y — u + — - logu J (3.35) 

. ( cosh 9 sinhG \ 1 i (uj + m\ 

^ = i smhe coshe j ' ® = 4 ^°H^J ^'-''^ 

(3.37) 



u 

and a constant C > 0. 



Proof. We write (3.7) symbohcally as 

X' = V X 

with a matrix potential V{u). According to (3.30) and the hypothesis |a;| > 
m, the eigenvalues of V are, for sufficiently large u, purely imaginary. More 
precisely, there is a transformation matrix B{u) with 

B-^VB = -iQa^ (3.38) 

and a suitable function Cl{u). Since the matrix potential V converges for 
u oo and has a regular expansion in powers of 1/u, we can choose B such 
that 

\B{u)\ < Co , \B'{u)\ < ^ (3.39) 

with a constant cq > 0. The transformed wave function {B~^X) satisfies the 
equation 

(B-^X) = \-in(u) - B-^ B'] (B-^X) . (3.40) 

au 

Hence employing the ansatz 

X = b(^ ^^7-(n? ) ^'^"^ = ^^"^ ^^-^^^ 

and using the bound (3.39), we obtain the inequality 

d 



du ^ 



< I I/I ■ (3-42) 



A short calculation shows that <I> has the explicit form (3.35), and that 
B{u) = A + 0{^) with A according to (3.36). The term of order 0{^) can 
be absorbed into i?no- 



F. Finster, N. Kamran, J. Smaller, S.-T. Yau 



43 



The inequality (3.42) can be used similar to (3.14) in Lemma 3.1 Namely, 
dividing by |/| and integrating yields for sufficiently large u the bounds 

J < I/Ml < L. (3.43) 

After substituting the upper bound for |/| into (3.42), one sees that /' is 
integrable. Thus / has a finite and, according to (3.43), non-zero limit, 

/oo := lim f{u) + . 

M— >0O 

Finally, the l/?i-decay (3.37) follows by integrating (3.42) backwards from 
■u = oo and employing the resulting bound in the ansatz (3.41). ■ 

In analogy to potential wall problems for Schrodinger operators, we call the 
function /oo in (3.33) corresponding to our fundamental solutions '^^^ the 
transmission coefficients, and denote them by j^^fi- 

Theorem 3.6. For every * G C^(IR x S'^f , 

(3.44) 

where the coefficients tab for |a;| < m given by 

tab = ^b,i ■ (3.45) 
For \lo\ > m, the tab o-fe given by the integrals 

where the functions ta are related to the transmission coefficients by 

= /+2e-^"-/oo2e'", *2(a) = -/+^e— + /^ie-. (3.47) 

The integral and the series in (3.44) converge in norm in the Hilbert space 

n. 

Proof. Our strategy is as follows. Choosing U2 so large that supp * C 
(—00,^2), Proposition 3.4 yields for t = the "completeness relation" 

1 

*(^) = ^ E / <*«r I ■ 

fcez -^-^ nez 
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This formula remains true when U2 is further increased, 

1 r°° 

^(^) = ZY. KT+ri^) <^«^^r I ^> , ^ > 0. (3.48) 



fcGZ """" nGZ 



Hence we can take the average over r in the finite interval [0, T] with T > 
and obtain, using Fubini's theorem. 



1 r°° r 1 /"^ 

feez -^-^^ nez 



(3.49) 



We shall first prove that the square bracket in (3.49) has a finite limit as 
T oo. Then wc will show that for T — > oo, wc can in (3.49) take the limit 
inside the integral and the series in (3.49). This will give a decomposition of 
the identity in terms of eigensolutions of H, from which the representation of 
the propagator (3.44) will follow immediately by inserting the phase factors 
ex.p{—iLot). 

Let us analyze the square bracket in (3.49). We can write '^u2+t ^ ^ 
linear combination of the fundamental solutions 

2 
a=l 

where the coefficients Cij^ must be chosen such that our Dirichlet-type bound- 
ary conditions are satisfied at U2 + t. Then the square bracket becomes 



T 
with 



1 ^ 

- / dr ^t^+r <<2+r I ^> = E *«^(^) <*^" I (3-51) 

^ -^0 a,b=l 



tab{T) = ^ [ Ca{T) Cb{T) dr . (3.52) 

^ Jo 

In the case \uj\ < rn, vj/^'^" and ^2'^" arc for large u exponentially decaying 
and increasing, respectively. Hence in order to fulfill the boundary condi- 
tions at u = U2 + T, the quotient C2(r)/ci(r) must go exponentially to zero. 
Moreover, our normalization conditions (3.17) and (3.31) imply that |ci(r)|^ 
must tend to one. We conclude that there is a constant ci with 



\ca{r) - Sa,l\ < Ci e' 

and so (3.52) converges for T — > 00 to (3.45). In the case > m, the funda- 
mental solutions are oscillating for large u, as described by Lemma 3.5. The 
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boundary conditions at U2, (2.21), imply that the following scalar product 
must vanish, 

where /ly^ are the transmission coefficients. Moreover, the normalization and 
phase conditions (3.17) and (3.32) yield that 

|cip + |c2p = 1. (3.54) 

The general solution to (3.53) is 

with a complex parameter D, which can be chosen so as to satisfy the nor- 
malization condition (3.54). We now substitute (3.55) into (3.52) and take 
$ as the integration parameter, 

U(T) = - / Cam cm ^ ■ (3.56) 

Using (3.35), one sees that (3.56) converges for T ^ oo to the average over 
one period, giving (3.46) and (3.47). We conclude that the bracket in (3.49) 
converges pointwise as T — *■ oo. 

Next we shall prove that in (3.49) we may take the limit T ^ oo inside 
the series and the integral, and that for the resulting limit the series and 
the integral converge in norm. The sum over k in (3.49) gives the decom- 
position into the eigenspaces of the angular operator id^, which we denote 
by Ti^. We may consider the situation on each such cigcnspacc separately, 
and thus assume that ^' G TC^. For the integral and the n-summation in 
(3.49) the situation is more difficult because the spectral decomposition of 
the Hamiltonian depends on U2, and because the eigenvalues A„ of A and 
corresponding eigenspaces depend on to. We first apply to (3.48) the opera- 
tor product A'^^H'^'^ with p,q > (with A as in (2.4), where the t-derivative 
in A is carried out by applying the operator —iH) and take the inner product 
with This gives 

/oo ^ 
— co^i ^'ni^) \<K'^+T I ^>l^ • (3.57) 

If we consider on the right side of (3.57) instead of |<^'^^"^|^'>p a mixed 
product with * G C^((— oo, U2) x S'^)^ CiH'' , we can in the integrand use 
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the inequality xy < ^(x^ + y^), and then apply (3.57) to obtain 



nez 



2q 



< 2 I -^^^ ^^'^ + <^ I i/^*^ *>) , (3.58) 
and this bound holds for all r > 0. 

Let e > 0. Then for any loq > 0, (3.58) yields for p = and q = I that 

'^'^ <$ I ><*'='^" I 

^ I ^U2+T^^ U2+T I ^ 



1 roo ^ 



^ ^{\mf + \mf) 



We choose so large that 

dw 



< 



(3.59) 



for all T > 0. This inequality allows us to restrict attention to uj in the finite 
interval [— u;o,u;o]. Next for a constant no > 0, we consider (3.58) for p = 1 
and q = 0. This gives the inequality 



^ r V 

|n|>no-^-'^o 



< - (<$ I ^2$> + <* I ^2*>) sup K {^) ■ 

^ a;e[-a;o,a)o], |n|>no 

Clearly A^(a;) — > oo for n — > ±00 uniformly in a; G [— a;o,wo]) and thus we 
can by choosing no sufficiently large arrange that 



E 



^0 duj 



TT 



<<B I ><\I/'^''^" I \I/> 



< 



(3.60) 



|n|>no''~'^° 

Putting together the estimates (3.59) and (3.60), we conclude that 



/oo 
- E 
-00 I I. 



E 



\n\>no 



LOO > ^ 



I ^ U2+T ^ ^ U2+T I 



< e (3.61) 
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for all r > 0. 



For a; in a finite interval and n in a finite set, we can take the average 
over T G [0,T] and take the limit T ^ oo using Lebesgue's dominated 
convergence theorem (notice that according to Lemma 3.1, ^^J^^ is on the 
support of $ and ^ uniformly bounded in r). This gives 



lim 

T^oo 



no 

E 

i=-no 
no 



-uiQ ^ Jo 



= J2 Yl *ar<* I *a'^"> <*&'^" I *> (3-62) 

n=-no-^-^o a,b=l 

with tab according to (3.45) and (3.46). Since e in (3.61) can be chosen 
arbitrarily small and uq oo, (Jq ^ oo as e \ 0, we obtain that (3.62) is 
true also for no = oo = wq, with absolute convergence of the integral and 
the series. Since $ can be chosen arbitrarily, we conclude that 

poo J 2 

= lim (3.49) = V / ^a'^" <*6'^" I *> ■ (3-63) 



— l^t^b 

a,b=l 



The estimate (3.58) for p = = q remains true if we take the average over 
T G [0, oo), and a homogeneity argument (as in the proof of the Schwarz 
inequality in Hilbert spaces) yields that 



nez 



-oo ^ 



J2 tab I *^'^"><*^'^" I 



a,b=l 



< ll$ll 11*11 



This bound shows that the integral and series in (3.63) converge in norm, 
and that * need not be an eigenvector of id^. We finally apply the unitary 
operator exp(— ziif) on both sides of (3.63) to obtain (3.44). ■ 



Notice that the coefficients tab given by (3.46) are bounded, 

\tab\ < - for \uj\ > m. 



1 

< - 
- 2 



(3.64) 



In the asymptotic region u — > — oo, (3.44) goes over to a Fourier represen- 
tation in terms of the plane- wave solution (3.8). A careful analysis of this 
limiting case gives additional information on the coefficients tab, namely 

til = - = ^22 for > m. (3.65) 

However, the non-diagonal elements ti2 and t2i remain undetermined. We 
shall not derive the relations (3.65) here, and will not use them in what 
follows. 
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4 The Decay Estimates 

Using the integral representation for the propagator of the previous section, 
we can now prove the decay of the probabihties. 

Proof of Theorem 1.1. According to the hypotheses of the theorem, the 
initial data ^'o is in L^((ri,oo) x S^,diJ,)'^. Since the transformation of the 
spinors (2.1) is smooth and involves factors A* and -s/r, we obtain for the 
transformed initial data *o that 

r-5 A-2 G L^((ri,oo) x . 

Equivalently, computing the volume element on the hypersurface t = const, 

r^ . 

— |*oP e L^{{ri,oo)x S^,drdcos-&dip)'^ 

Transforming to the variable u, (2.10), one sees that 'I'o is in the Hilbert 
space with scalar product (2.16), £ For simplicity, we again omit the 
hat in what follows. 

Let e > 0. Since the wave functions with compact support are dense in 
Ti, then for this e, there is a ^'7 G C°°((ri,oo) x S'^)^ such that 

W'Si - ^'oll < e . (4.1) 

For the Cauchy problem with initial data we have the integral represen- 
tation of Theorem 3.6. Since the series in (3.44) converge in norm, we can 
choose ko and no such that 

||*iko,no-*/|| < e, (4.2) 

where ^'jto,no is defined by 

*feo,no(^) = - E E/ E*'r*^"(^)<*^"i*/>- (4.3) 

fc=-feo n=-no a,b=l 

Consider the integrand in (4.3) for fixed A; and n, 

2 

E *a'^''(a;) <*6'^'' I *7> • (4.4) 

a,b=l 

iFrom (3.64) and the estimates of Lemma 3.1, one sees that (4.4) is bounded, 
locally uniformly in x and to. Thus the norm convergence established in 
Theorem 3.6 implies that (4.4) is in L^(IR, (D^) as a function of u, with 
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an L -bound locally uniform in x. Hence its Fourier transform is L°° in t, 
locally uniformly in x. Furthermore, the Riemann-Lebesgue lemma [8] yields 
that its Fourier transform tends to zero as t — > oo, pointwise in x. Since 
(4.3) involves only finitely many terms, we conclude that the solution of the 
Cauchy problem with initial data '^ko,noj 

1 *:o "0 poo 2 
^ k=-ko n=-no a,b=l 

(4.5) 

is L°° in t locally uniformly in x, and lim(_^.oo '^no,ko(^^ x) = for all x. 

Choose Ks^R as in the statement of the theorem. Since the metric and 
Dirac matrices in the probability integral (1.3) are smooth and bounded 
on the compact set K^^ji, the corresponding bilinear form is continuous on 
H, i.e. there is a constant c depending only on S and R such that for all 
*i,*2 eH, 

/ (¥ij^^2) i^j dfj, < c ||*i|| ||*2|| ■ (4.6) 

JKs,R 

The solution to our original Cauchy problem is obtained by applying the 
unitary operator e:x.p{—itH) to ^'o, 

= e-^*^ ^ko,no + e-'*^ i^i - ^ko,no) + e~**^ (*o - */) 

= *fco,n«(*) + e-**^ i^I - ^ka,no) + e'"'" (*o - */) , 

where ^t'fccnol^) has the integral representation (4.5). We substitute this 
formula for '^{t) into the probability integral, multiply out, and apply the 
estimate (4.6) as well as the unitarity of exp(— ■itii") together with (4.1) 
and (4.2). This gives the inequality 

/ {^'j^^){t,x) Vj diJ, 

< / (*feo,no7^*fco,no)(*> x) Uj dfi + 4c^e^ + Acs W^W . 

We showed above that the integrand in the last integral is uniformly bounded 
and tends to zero pointwise as t — oo. Thus the integral converges to zero 
according to Lebesgue's dominated convergence theorem. ■ 

We remark that in the spherically symmetric case, the analytical method 
given above to prove that ^ 7i \s an alternative to the nice geometric 
argument by Kay and Wald [7] , who use the causal propagation property and 
a discrete symmetry of the maximally extended space-time at the bifurcation 
2-sphere. 
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A Nondegeneracy and Regularity of the Angular 
Eigenfunctions 

In this appendix, wc shall consider the angular equations (2. 6), (2. 8). As 
explained in [1, Appendix A], it is useful to write (2.6) as an eigenvalue 
equation in A, 

. , . / —am cost? C- \ /. -.x 

AY = XY with A = { ^ a . (A.l) 

y —C+ am cosv J ^ ' 

Proposition A.l. For given k and X G (^{A), there is at most one eigen- 
solution of (A.l), which we denote by Y^ , i.e. 

AY^ = XY'' . (A.2) 

By continuously varying the parameter u), the eigenvalue equation (A.2) can 
he extended to all values of lo E IR. Both X and depend smoothly on uj. 

Proof. The two fundamental solutions of (A.2) behave near ■!? = like 

= + o{^''), 0(1?^=)) and = (0(7?-^-^), 1}-''-^ + o{^-''-^)) , 

respectively. Depending on whether A; is > or negative, the second or first 
fundamental solution diverges in the limit 'd ^ 0. In [1, Appendix A] it was 
shown that the eigenfunctions are bounded on S"^ and smooth except 
at the poles. Thus we can rule out one of the fundamental solutions and 
conclude that (A.2) has at most one solution. 

Note that the solutions of (A.2) are the eigenvectors of A restricted to 
the eigenspace of the operator id^p with eigenvalue k, which we denote by H''. 
Since the terms involving iv in (A.l) arc a relatively compact perturbation, 
standard perturbation theory [6] yields that the spectrum of A\-}^k depends 
continuously on a;. As no degeneracies occur, each eigenvalue A gives rise 
to a unique continuous family of eigenvalues A(a;). Standard perturbation 
theory without degeneracies [6] then yields that A(a;) and the corresponding 
eigenvector Y'^{u}) depend smoothly on 00. ■ 
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